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Abstract 

We continue our studies of spherically symmetric self-similar solutions 
in the SU{2) sigma model coupled to gravity. For some values of the 
coupling constant we present numerical evidence for the chaotic solution 
and the fractal threshold behavior. We explain this phenomenon in terms 
of horseshoe-like dynamics and heteroclinic intersections. 

1 Introduction 

This is the fourth paper in a series aimed at understanding of the structure of 
the self-similar spherically symmetric wave maps coupled to gravity. In the first 
two papers P] G] it was shown that for the coupling constant a < | there exists 
a countable family of solutions that are analytic below the past light cone of 
the central singularity. It was also shown that there are two generic types of 
solutions beyond the past light cone. The threshold between these two generic 
types resembles many aspects of critical phenomenon and has been studied in 
P] for the coupling constant a < 0.426. It was found that for this coupling the 
threshold solutions are the periodic solutions. Now, we extend studies of the 
threshold solutions for larger coupling constant and examine the structure of 
the bifurcations in this parameter. This article is mathematically self-contained 
but for physical background we refer the reader to 0] . 

2 Setting the problem 

Under the assumption of spherical symmetry and self-similarity the Einstein's 
equations with wave map matter reduce to a system of autonomous equations 
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W = -l + a(l ~W 2 )F' 2 , (1) 

A' = -2aAWF' 2 , (2) 

, , , w sin(2F) , 

(AF ) = ^i-i, (3) 

subject to the constraint 

1 - A- 2a sin 2 (F) + aAF' 2 (W 2 - 1) = 0. (4) 

Physically the functions W and A parameterize the metric and the function F 
parameterizes the corotational wave map. a stands for a dimensionless coupling 
constant (a = corresponds to no gravity). We are interested in solutions of 
I^OJ starting at xq (we choose xq = 0) with the following initial conditions 

W(0) = 1, 4(0) = 1- 2a, F(0) = |, F'(0) = b (5) 

where b is a free parameter (since the system has reflection symmetry F — > — -F 
we can take b > without loss of generality). Solving Q for A and using 
a new variable D = F' we rewrite equations ^BJl as the three dimensional, 
autonomous, first order system (a < |) 

W = -1 + a(l - W 2 )D 2 , (6) 

£>' = 2^ + Bin(2F) / (7) 
-l + 2asin(F) 2 \ l-W 2 )' W 

F' = D. (8) 

with the initial conditions 

W(0) = 1, D(0) = b, F(0) = f . (9) 

In the following we will refer to solutions of ((HHSl satisfying the initial conditions 
© as b — orbits. It was shown in [2] that for a S [0, \ ) b — orbits exist locally 
and are analytic in b and x. 



3 Types of solutions 

It was shown in that for small values of b the b — orbits exist up to the 
future light cone (W = —1) of central singularity, while for large values of b the 
b — orbits develop a sort of singular apparent horizon beyond which they cannot 
be continued (W = +1). These two types of solutions, call them type A and 
type B solutions, respectively, can be shown to form open sets [Hj- Hence, there 
must exists solutions which are not type A or B, call them type C solutions. 
For a given a there exists a critical value b* (a) such that solutions with b < b* 
are of type A, solutions with b > b* are of type B, and the solutions with b = b* 
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is of type C; in other words, we have a bistable behavior with two generic final 
states A (W = -1) and B (W = +1), and the separatrix C (\W\ < 1 for all 
x). We use the same terminology to distinguish attractors A (W = — 1) and B 
(W = +l). 

It was shown in pj] that type C solutions are asymptotically periodic for 
a < 0.426. We will show in this paper that this remains valid for a > 0.482. 
However, we will also show that for a € (0.426, 0.482) (approximately) the 
transition between types A and B does not occur at a single point but in a 
narrow interval {bl ninl b^ nax ) and the behavior of type C solutions is chaotic. 
This chaotic behavior is the main point of our interest. 

4 Type C solutions 

Numerical data show that for general initial conditions the whole phase space 
contains two basin sets (of attractors A and B) separated by a two dimensional 
surface. Type C solutions lay in this critical surface. In other words all initial 
points in the critical surface give type C solutions. The critical surface is an 
invariant submanifold. 

We are especially interested in asymptotic structure of type C solutions. 
Hereafter we will skip initial part of C solutions to simplify all figures. Let 
us call our asymptotic structure C-attractor. In fact, it is the attractor of 
codimension one (two stable and one unstable directions). The structure of 
C-attractor for a < 0.426 has been found in p] it is the periodic solution. In 
other words, b — orbits starting close to the critical b* follow periodic C-attractor 
for some time and then run towards A or B attractor (see figure El . 

5 The bifurcation diagram 

Let us consider the initial curve in the form For each a, using the bisection 
procedure, we can find a critical values of b* parameter which determine inter- 
section of the initial curve with critical surface (figure |2J . It has been noted in 
PQ that for large values of a ~ 0.43, b* (a) is not a single point but there is a 
structure hidden in a narrow interval (b* m j„ (a) , b* max (a)). The solutions here 
are sensitive to initial conditions infinitesimally small change of b parameter can 
result in a completely different final state. We have confirmed these numerical 
results. Since the magnitude of interval (6* m , n (a), b* m ax{ot)) is of order 10~ 9 
we used quadruple machine precision. In addition, the initial conditions are 
given at a singular point. We are interested only in asymptotic structure of type 
C solutions so it is more convenient to reformulate initial curve to start directly 
from C-attractor 

W(0) = 0, D(0) = c, F(0) = 0, (10) 

where c is a function of b such that \c* max - c* min \ » \b* max - b* mm \. The 
function c(a) is qualitatively similar to b(a) but now at a approximately equal 
to 0.426 we can resolve better the new structure in the diagram (figure EJ|. This 
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is a bifurcation point. Numerical results suggest that this new structure is self- 
similar (figures IH1 EI HI so for each a we have infinitely many critical points. 
This suggests that the basin boundary is fractal. 

There is not enough numerical data to calculate the fractal dimension from 
the definition but the phenomenon of fractal basin boundaries is well known so 
we have another tools to do this. We use the procedure suggested in which 
goes as follows. 

We call a point e-uncertain if perturbations with e <C 1 in opposite directions 
give points which belong to different basin sets. It was shown in that the 
number of uncertain points in the phase space / scales as 



/(e) ~ e a . (11) 



Therefore, in the limit we have 



lim = a. (12) 

e^o me v ' 

It was shown in that the uncertainty exponent a is related to the capacity 
dimension d in D dimensional phase space by the formula 

d=D-a. (13) 

For a — 0.4264 we estimate numerically the capacity dimension to be 
d = 0.337 ± 0.003 (see figure^. Since C-attractor lays in critical surface with 
infinite area one should expect complicated dynamic. Figure shows only one 
dimensional relation. It would be interesting to look at the whole basin bound- 
ary. We can do this by taking slices of the phase space (see figure EJ . 



6 C-attractor 

We found C-attractor in this interval of a numerically (figures El and EJ . The 
problem is not trivial because C-attractor lays in the unstable manifold (two sta- 
ble and one unstable directions) . We use straddle-orbit method due to Battelino 
et al. [f3 and H. Nusse and J. Yorke 0. 

This procedure allows us to pursue the unstable orbit laying in the basin 
boundary in principle forever. It can be viewed as a series of bisections and 
goes as follows. At initial time we choose two points Pa(x = 0) and Pb(x = 
0) which lead to different attractors A and B and perform bisection until the 
distance between the iterates -Pa(O) and -Pb(O) is less a prescribed S. Next we 
integrate the equations numerically starting from the current -Pa(O) and -Pb(O) 
until the distance between the trajectories exceeds 5. When this happens at 
some time x we stop the integration, assign the points Pa(x) and Pb{x) as 
current representatives and repeat the bisection. Iterating this procedure one 
can progressively construct a trajectory staying within a distance S from the 
codimension one stable manifold. 
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Using straddle-orbit method we found two loops and a sequence of the wind- 
ings f figure Hot . It is in a way similar to the Lorenz attractor but the nature of 
C-attractor is different. The standard methods like Lyapunov exponents, test 
of the subsequent maxima, Poincare sections of the attractor do not give good 
insight into the problem although behavior of maximal Lyapunov exponent in- 
dicate the chaotic region (see figure 7 in pjj). The best evidence of chaotic 
properties provides the autocorrelation function which exhibits a typical beha- 
vior (figure EJ- It shows that we are able to predict short range behavior of 
the system (inside single loops) with good accuracy but long, deterministic se- 
quence of windings is uncorrelated. We can show this numerically up to almost 
any finite period with help of the symbolic dynamics and information theory. 
Symbolization 1 allows us to check periodicity of our solutions, calculate Shan- 
non entropy or EMC (Effective Measure of Complexity given by Grassberger 
|B]) and show sensitivity to initial conditions. 

7 The horseshoe-like dynamics 

The results presented in previous subsections are based on numerical calcula- 
tions. It is well known that chaos and fractal structures could appear in dynami- 
cal systems due to the discretization errors [9|. In order to make confident that 
our result is not a numerical artifact we checked (in the bifurcation diagram) 
the dependence of the size of chaotic windows on machine precision. We have 
performed this calculation with positive result: the size of the chaotic windows 
does not tend to zero with increasing machine precision (figure El- Unfortu- 
nately there is not enough computational efficiency to use more detailed lattice 
and estimate the limit. 

In summary, we have strong numerical evidence for chaos but we need a 
qualitative argument. Here enters the horseshoe. 

Now, we will show that there is two dimensional section of the basin bound- 
ary which has a fractal structure. 

Let us consider a two dimensional one to one and continuous map M. We 
choose two sets A and B such that A C A and B C B where A, B are two basin 
sets. Of course we have Ad B — 0. We choose a test set I (gray rectangle in 
figure n"8ll and iterate M(I). If M(I) and / cross in the way presented in figure 
EH we could say something more about structure of set / (see figure EJ. There 
are subsets in I which go to A or B. By D we denote this part of domain for 
which we are not sure where it goes to. The repetition of this procedure will 
reveal more and more complicated structure of the test set. In the limit we will 
get a structure similar to a fat Cantor set. 

With the help of the Poincare map we could reduce three dimensional flow 
to two dimensional map (Poincare map is one to one and continuous as map 
M). Test set and its intersection satisfying described conditions are presented 
in figure IT51 

*It will be described elsewhere. 
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It follows from the foregoing that there is a surface in the phase space whose 
intersection with the basin boundary has a fractal structure. This, together 
with numerical data implies existence of the chaotic solution. We believe the 
mathematically rigorous proof can be done by applying the method presented 
in [Ej. 

8 Route to chaos 

There is also a parallel description of the phenomenon presented in the previous 
section which helps to follow "the route to chaos" in Poincare picture. 

The system ijliHEJl has reflection symmetry F — > — F so for each initial value 
of c* also — c* gives critical solution. We take a slice of the whole phase space 
and we denote by N the map defined by subsequent intersections of the flow with 
the chosen slice. This corresponds to four period two points for small positive a 
(intersection of two loops with the surface) . Map N 2 reduces four periodic points 
to two saddles and gives us proper Poincare section. These saddles lay in the 
basin boundary. The intersection of the basin boundary with Poincare surface 
gives stable manifolds of the saddles for N 2 . If we choose Poincare section with 
F = const, we can use existing symmetries to determine unstable manifolds. 
The derivatives W , D' are invariants of the reflection W — ► —W, D — > —D. 
It follows from this that stable manifolds under this transformation correspond 
to unstable manifolds (see figure I Kill . The crossings determine positions of the 
saddles. The shapes of the manifolds and positions of the saddles depends on 
the parameter a. 

For small a, b — orbits have been found analytically (see jHj)-) We know from 
this three dimensional analysis that for a — > periodic solutions run towards 
infinity (as we can see in figure^ and their period shrinks to zero. That means 
in terms of the Poincare picture that the distance between saddles is unbounded 
for a — > 0. 

If we go in the opposite direction (growing a) the distance between the 
unstable manifold of the first saddle and the stable manifold of the second 
saddle decreases. At the bifurcation point they collide (a ~ 0.426) and the 
heteroclinic intersection arises (see figure H?)jl . This implies an infinite number 
of intersections and a chaotic nature of C-attractor. For larger values of a 
(a ~ 0.482) the heteroclinic intersection disappears and C-attractor becomes 
periodic again. 

The bifurcation of the basin boundary seems to be different from the one 
"at point" described for the Henon map. We have observed that manifolds 
are deformed just before bifurcation point (a = 0.4244, figure Il7ll . Further 
investigation is needed to clarify this fact. 







9 Final remarks 



The main result of this paper is the qualitative and quantitative argument for 
the existence of a chaotic solutions of the equations 111131 . We have skipped 
the physical background but we point out that these equations are reduced 
Einstein's equations. In addition, chaotic solutions separate two generic types 
so in this sense they are critical. As far as we know this is the first example of 
fractal critical behavior in the context of Einstein's equations jllj . 
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Figure 1: The function W - the type C solution (periodic). The type A solution 
(c = 2.36134), the type B solution (c = 2.36135). All calculations for a = 0.38. 
The type A, B solutions follow periodic C-attractor for some time and then run 
towards A and B attractor. 
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Figure 2: The bifurcation diagram, lattice Aq = 10~ 3 i Afe = 2 ■ 10~ 3 . 
The critical curve b*(a) separates type A solutions from type B solutions. 
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Figure 3: The bifurcation diagram, lattice Aq = 10 -4 i Ac = 2 ■ 10~ 4 . 
The critical curve c(a) separates type A solutions from type B solutions. 
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Figure 4: The bifurcation diagram (figure EJl enlarged 1000 times , lattice 
Aa = 10~ 4 and Ac = 2 ■ 10 -4 . The complex structure of the basin boundary is 
visible. 
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Figure 5: The bifurcation diagram (figure EJ) enlarged 200000 times, lattice 
Aa = 10~ 5 i Ab = 10~ 6 . The complex structure of the basin boundary is 
visible. 



13 



-9 -8 -7 -6 -5 

ln(E) 



-4 -3 



Figure 6: The number of uncertain points in the phase space (for a — 0.4264) 
as a function of the perturbation parameter e (logarithmic scale). The linear 
fit gives the slope a = 0.663 ± 0.003 from which using 11211 and 111 -it we get the 
capacity dimension d = 0.337 ± 0.003. 



14 



F=0 




-0.3 0.3 0.6 
W 

F=0.87 




W 



Figure 7: The basin boundary for a 
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0.43. The complex structure is visible. 
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Figure 8: C-attractor in phase space (W, D, F) (for F we identify the intervals 
[-§ + irk, f + irk] where k € N), a = 0.43. 
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Figure 10: The projection of C-attractor, a = 0.43 (for F we identify the inter- 
vals [— f + 7rfc, ^ + irk] where k £ N). Some similarities to the Lorenz attractor 
are visible. 
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Figure 11: The correlation function; the comparison of the periodic solution 
a = 0.42 with chaotic a = 0.43 (only short time correlation exist). 
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Figure 12: The size of chaotic window as a function of the numerical precision. 
Since the result does not tend to zero with increasing machine precision we 
believe that the observed chaotic behavior is not a numerical artifact. 



20 




Figure 13: Sets: A, B, test set I and M(I) (S-shape). We have: A C A and 
B C B where A, B are two basin sets. The horseshoe in the crossing M(I) and 
/ signals fractal boundary between A and B. 
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Figure 14: The structure of the set / determined after 0, 1, 2 iterations of the 
map M. The repetitions of this procedure converge to a part of the boundary 
(between basin sets A and B) and will lead to a fat Cantor set in the limit. 
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Figure 15: The intersection of the test set I (bounded by straight lines) and 
its first image by Poincare map. We have plotted in the background the basin 
boundary (visible as four, almost parallel curves). These curves correspond to 
the stable manifolds of the periodic points. The first iteration of the test set I 
corresponds to the unstable manifolds. We observe horseshoe as in figure IT7T1 
This is a "first approximation" of the heteroclinic intersection (figure IT7)j) . 
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Figure 16: Two saddles and the creation of the heteroclinic intersections of 
stable and unstable manifolds. 
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Figure 17: The bifurcation of the basin boundary, F=1.0. The strong deforma- 
tion is observed (a = 0.4244) before a chaotic behavior appears (a ~ 0.426). 
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